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admit propagating singular solutions are introduced and discussed. 
These singular solutions result of the competition be- 

£h ! tween nonlinear and nonlocal processes in various familiar vector 

spaces. Several examples are given. The motivating example is 
the directed self assembly of a large number of particles for tech- 
q : nological purposes such as nano-science processes, in which the 

■ particle interactions are anisotropic. This application leads to the 

derivation and analysis of gradient flow equations on Lie algebra 
valued densities. The Riemannian structure of these gradient flow 
equations is also discussed. 
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1. Introduction 

Many physical processes may be understood as aggregation of individual 
'components' at a variety of scales into a final 'product'. Diverse examples of 
such processes include the formation of stars, galaxies and solar systems at 
large scales, organization of insects and organisms into colonies at mesoscales 
and self-assembly of proteins, nanotubes or micro/nanodevices at micro- and 
nanoscales [I]. Some of these processes, such as nano-scale self-assembly of 
molecules and particles are of great technological interest 0, Ej. In most 
practical cases, the assembling pieces do not have spherical symmetry, and 
the interaction between those pieces is not central, i.e., it is dependent not 
only on the density of the particles, but also on their mutual orientation. A 
recent example, motivating this paper, is self-assembly of non-circular float- 
ing particles (squares, hexagons etc.) [H O, Ej. Due to the large number of 
particles involved in self-assembly for technological purposes, (10 9 - 10 12 ), 
the development of continuum descriptions for aggregation or self-assembly 
is a natural approach toward theoretical understanding and modeling. 

Progress has been made recently in the derivation of continuum evolutionary 
equations for self-assembly of molecules evolving under a long-range central 
interaction potential. The resulting continuum models emerge in a class of 
partial differential equations (PDEs) that are both nonlinear and nonlocal 
[H [S] • Nonlinearity and nolocality may combine in an evolutionary process 
to form coherent singular structures that propagate and interact amongst 
themselves. Thus, self-assembly is modeled in the continuum description as 
the emergent formation of singular solutions. Classic examples of equations 
admitting such singular solution behavior include the Poisson-Smoluchowsi 
(PS) equation 0, Debye-Hiickel equations PH] and the Keller-Segel (KS) 
system of equations [TT]. PS and KS are gradient flows for densities that 
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collapse to subspaces in finite time. Their respective collapses model the 
formation of clumps from randomly walking sticky particles (PS) and the 
development of concentrated biological patterns via chemotaxis (KS). For re- 
cent thorough review and references, see [12]. A variant of these equations 
of potential use in modeling formation of clustering in biology (e.g., forma- 
tion of herds of antelope, or schools of fish) was introduced in Bertozzi and 
Topaz [13]. Dynamics and collapse of self-gravitating systems has been stud- 
ied by Chavanis, Rosier and Sire [H]. Finally, aggregation of particles with 
the applications to directed self-assembly in nanoscience was considered by 
Holm and Putkaradze (HP) in (TJ, [SJ . The singular solutions of the latter HP 
equation emerge as delta functions supported on a subspace from smooth 
initial conditions. These clumps or "clumpons" then propagate and interact 
by aggregating into larger (more massive) clumps when they collide. See 
also [T5] for a thoughtful discussion of singularity formation in fundamental 
mathematical terms. 



Remarkably, the aggregation and emergence of singular solutions from smooth 
initial data due to nonlinearity and nonlocality need not be dissipative. A 
Hamiltonian example is provided by the Camassa-Holm (CH) equation for 
shallow water waves in its limit of zero linear dispersion [TB]. In this limit, the 
CH equation describes geodesic motion on the Lie group of smooth invertible 
maps with smooth inverses (the diffeomorphisms, or diffeos, for short) with 
respect to the H 1 Sobolev metric. Being an integrable Hamiltonian system, 
the CH equation has soliton solutions which emerge in its initial value prob- 
lem. In its dispersionless limit these CH soliton solutions develop a sharp peak 
at which the derivative is discontinuous, so the second derivative is singular. 
At the positions of these propagating peaks (the peakons), the momentum 
density is concentrated into delta functions. The CH peakons propagate and 
interact elastically by exchanging momentum; so they bounce off each other 
when they collide. 

CH arises as an Euler-Poincare (EP) equation from Hamilton's principle de- 
fined using a Lagrangian which is right-invariant under the action of the 
diffeos on their own Lie algebra (the tangent space at the identity) [Hj. In 
the EP approach, variations in Hamilton's principle of the vectors in linear 
representation spaces of the group of diffeos are induced from variations of 
the diffeos themselves. Hence these variations in the EP Hamilton's principle 
arise as geometrical operations, primarily as Lie derivatives of elements of 
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the appropriate linear representation spaces. Moreover, evolution by CH - 
or more generally by any of the class of EP equations on the diffeos (EPDiff) 
- is by coadjoint action of the diffeos on the momentum density. This in- 
finitesimal action is defined as the Lie derivative of momentum density by 
its corresponding velocity vector field. The momentum density for a singular 
solution of CH is supported on its singular set, regarded as an embedded 
submanifold [TBI • 



CH and other equations in the class of EPDiff equations hint that the actions 
of diffeos may be taken as a paradigm for the development and propagation 
of singularities in other geometrical quantities, or at least would provide a 
framework for studying the process and discovering additional examples of it. 
In fact, the HP equation was discovered by following Otto [12] in formulating 
the process of nonlinear nonlocal evolution using a variational principle in 
which the variations were induced by the infinitesimal actions of diffeos on 
densities f7\ IB]. Just as in Hamilton's principle for the EP equations, this 
infinitesimal action is the Lie derivative with respect to a vector field and for 
HP the vector field was related by Darcy's Law to the flux of density. 

The derivation of evolution equations for density in the framework of Darcy's 
law is now well established: The velocity of a particle is taken to be propor- 
tional to the force acting on it, and the conservation law for density p readily 
establishes a PDE for its evolution. The Darcy law evolution equation for 
density corresponds to conservation of the n-form in n-dimensional space 
pd n x along characteristics of a velocity determined from the density and its 
gradient. Physically, the conservation law for pd n x means preservation of the 
number of particles, or mass, in an infinitesimal volume. 

The density evolution equation closes if the potential of interaction between 
the particles depends only on their relative position. This framework is sim- 
ple and attractive. However, the energy of some physical systems depends 
strongly on additional geometric quantities, such as the mutual orientations 
of their pieces. Examples of such systems are numerous and range from 
micro-biological applications (mutual attraction of cells, viruses or proteins), 
to electromagnetic media (dipoles in continuous media, orientation of do- 
mains), to interaction of living organisms [T]. 
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Summary of the paper. This paper provides a new geometric framework 
for continuum evolutionary models of particle systems whose interactions 
depend on orientation. Interestingly enough, this framework is completely 
general and can be applied to the evolution of any geometric quantity. The 
problems which involve both density and orientation are perhaps the most 
intricate, and for pedagogical reasons we shall treat them towards the end of 
the paper, after we have first illustrated the general framework by treating 
the familiar cases of scalars and 3-forms (particle densities). 

After this preparation establishing the pattern for the results in simpler cases, 
we shall derive the covariant evolution equation ( 19.71 ) for quantities which de- 
pend on both particle densities and orientations. In this framework, one 
may use the same principles in deriving evolution equations for any geomet- 
ric quantity. Of course, the forms of the various evolution equations will 
strongly depend on what geometric quantity is being considered. 

Once the equations of motions for the geometric quantities are derived, we 
shall concentrate on the establishment of singular analytical solutions of these 
equations. Earlier work for the case of density evolution [3, [H] demonstrated 
that when such solutions exist, they play a dominant role in the dynamics. 
Thus, we seek additional examples of the formation of coherent, propagating 
singularities in solutions of evolutionary PDEs. These singularities arise as 
a result of the competition between nonlinear and nonlocal processes in var- 
ious familiar vector spaces. We follow the same strategy as [3, [H] in taking 
variations of a free energy that are induced by infinitesimal actions of diffeos 
on these vector spaces. To close the equations, we also introduce a geometric 
analog of Darcy's Law, based on the dual representation of these actions. 
This approach yields several other nonlinear nonlocal PDEs whose solutions 
form coherent propagating singular structures in finite time. The present 
approach should be useful whenever geometry figures in the formation and 
evolution of moving singularities. However, the development of this approach 
is only beginning and the full extent of its applications remains to be proven. 
Besides the applications mentioned earlier, see for example, the spontaneous 
formation of singularities (self -focusing) in the spin wave turbulence |2"0] . 

Plan. The plan of the paper is the following. Section 2 reviews earlier work 
on the gradient flow evolution for density and formulates the geometric order 
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parameter (GOP) equation. This formulation summons the diamond oper- 
ator from differential geometry, whose properties are discussed in Section 3. 
Section 4 introduces a necessary condition for the existence of singular so- 
lutions of GOP equations in the general case. Applications require explicit 
calculation of the diamond operator for each geometric quantity, which is 
accomplished in Section 5. Section 6 derives explicit expressions for GOP 
equations of motion and Section 7 gives examples of singular solutions. Fi- 
nally, Section 8 is devoted to the equations of motion for orientation densities 
whose Lie-algebra-valued singular solutions are called gyrons. We also com- 
pare the results of our theory to recent experiments on self-organization of 
oriented particles. 



We consider continuum evolution of the macroscopic state of a system of 
many particles at time t and position x that is defined by an order parameter 
K(x,t), which take values in a vector space V. The vector space has a dual 
space V*, defined in terms of the L 2 pairing 



For example, scalar functions are dual to densities, one-forms are dual to 
two-forms and vector fields are dual to one-form densities. 

In addition, we assume that the physical situation dictates a free energy, 
which is a functional of the order parameter expressed as E[k], where square 
brackets [ • ] denote dependence which may be spatially nonlocal. That is, 
E[k] may also be a functional of k; for example, it could depend on the 
spatially averaged or filtered value of k defined later. Hence, the variation of 
total integrated energy is given by the pairing, 



where dot ( • ) denotes the appropriate pairing of vector and covector indices 
of k to produce a density, an n— form (denoted as A") which then may be 
integrated to yield a real number. Of course, n is the dimension of the space. 
In this paper, we shall concentrate on the cases n = 1,2,3. In this setting, 
we seek evolution equations for the order parameter K,(x : t) that 

(1) respect its vector space property k 6 V; 

(2) reduce to gradient flows when k is a density; and 

(3) possess solutions that may aggregate the order parameter k into "clumpons 



2. Problem statement 



(•,.): V x V* h-> M . 
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(quenched states) which propagate and interact as singular weak solutions. 

Thus we seek evolution equations whose solutions describe a geometric or- 
der parameter supported on embedded subspaces of the ambient space. For 
example, these solutions may be spatially distributed on curves in 2D, or 
on surfaces in 3D. In fact, we seek evolution equations for which these em- 
bedded singular solutions are attractors which emerge even from arbitrary 
smooth initial conditions, as in the case of the emergent singular densities 
for the nonlocal gradient flows studied in [7), EJ ■ Our ultimate goal is to find 
classes of equations that are relevant for the evolution of a macroscopic order 
parameter that may be of potential use in the design of directed self-assembly 
processes in nanoscience. 

Previous work focused on emergent singularities in an order parameter den- 
sity [7], [H] • These singularities correspond to the formation of dense clumps of 
maximum possible density in self-assembly of nano-particles. Two different 
cases were distinguished in the model. The first case arises when the mobility 
of the particles p always remains finite, in which case the evolution forces the 
particles to collapse into a set of (^-functions. These c^-functions can be under- 
stood as a large-scale view of individual clumps. The second case arises upon 
modeling the clumping process in greater detail, by limiting the mobility so 
that it will eventually decrease to zero at some maximal value of averaged 
density. This leads to formation of patches of constant density. (See also [21] 
for the case of density-dependent mobility.) In the next section, we shall cast 
the results of previous work into the geometric framework taken here in the 
derivation of the equations. 

2.1. An example: the HP equation for the gradient flow of a den- 
sity. 

Holm and Putkaradze (7J, E] derived the HP gradient flow, whose singular 
densities (the clumpons) emerge in finite time for any smooth initial condi- 
tions possessing a maximum density. The HP equation is the gradient flow 
of a density p G A" given by 



Equation ( 12.11 ) has a clear physical meaning as follows. A potential at a 
given point 5E/5p gives rise to the force V(SE/5p). The particles move with 
covariant velocity T] proportional to the force applied to them. The mobility p 
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gives the coefficient of proportionality. Thus, the instantaneous contravariant 
velocity at a given point is 

(2.2) rj= {pVSE/Spf. 

The familiar musical operations sharp (jj) and flat (b) raise and lower vector 
indices, respectively, thereby mapping a covector into a contravector and 
vice versa, as needed for the operations of divergence, Lie derivative, etc. to 
make proper mathematical sense. In particular, the contravariant velocity u 
advects the particle density according to 

(2.3) ^ = -£ { ^sE )iP , 



which provides the geometric meaning of the HP evolution equation ( 12.11 ) 



In the traditional framework, each extensive quantity is assumed to be con- 
served, so its density per unit volume p satisfies the continuity equation 

(2.4) | = - d - J . 

where J is the flux density of the conserved quantity. This flux density is 
assumed to be proportional to the density of the extensive quantity. It is 
also taken to be proportional to the gradient of the thermodynamic conju- 
gate variable 5E/5p in the First Law for specific energy E, modified to be 
contravariant so that its divergence may be properly defined. Thus the name, 
"gradient flow." Finally, a phenomenological characterization of the freedom 
of motion of the particles is introduced, called their "mobility" fi[p], which 
transforms as p. Thus, one finds the Darcy Law for flow through a porous 
medium, for example, as, 

SE\I 



(2.5) J = ^ V fo) 



Sp 

By this thermodynamic reasoning one recovers the HP equation ( 12.11 ). from 
the continuity equation ( 12.41 ) and the definition of mobility in Darcy's Law 
( 12.51 ). As a consequence, the energy E[p] evolves according to 

dE I ( _SE\l / _SE\\ 



" v ^)'(' v 



dt \ V Sp J V bp 

5E 



(2.6) = / pp[p] 



V 



d n x = I —-\J\ 2 d n x 



Sp J pp[p] 



FORMATION AND EVOLUTION OF SINGULARITIES IN ANISOTROPIC GEOMETRIC CONTINUA 



n 



Note that for attracting particles, E < 0, so the absolute value of energy 
decays in time d\E\/dt < when pp\p] > 0. 

Holm and Putkaradze 0, E] analyzed the solutions of the gradient flow (12.11) 
when the mobility density p[~p] depended on the spatially averaged density 
p = H * p, defined as the convolution of density p with the kernel H. The 
latter was chosen to be the Green's function for the Helmholtz operator. The 
energy density was chosen in [71 [8] so that 5E/5p = G * p, where G was a 
Helmholtz kernel whose scale length is larger than that for H. Thus, 5E/Sp 
depended on the average density at the lengthscale of G, rather than the 
pointwise density. 

With these choices, the solution p of the HP gradient flow, paired with a 
smooth scalar test function satisfies 



Under the the L 2 pairing (• , •), the third line matches with the derivative of 
a delta function; so the HP gradient flow admits singular solutions. 

2.2. Singular solutions for the HP gradient flow of a density. 

Provided p[~p] and 5E/5p are sufficiently smooth, the HP equation (12.11 ) ad- 
mits solutions p that are sums of delta functions supported on N manifolds 
of dimension K ai a = 1, 2, 3, . . . , N, with coordinates s embedded in R 3 . 
Namely, 




(2.8) 
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Integrating this solution ansatz for density p(x, t) against a smooth test func- 
tion (j) yields the dynamics of the parameters p a {t, s) (which are also densities) 
and the positions q a (t,s) G M 3 , as follows. After integrating by parts, one 
finds (suppressing obvious subscripts on s a and s&) 

r N r 

(<f>, pt + V- (p//(p)V(G*p))) = /#*)£ / Pa{t,s)S(x-q a {t,s))d 3 x 

J a=l 

» N . N . 

Pa (qq ~ ^ / Pb(t,s)ii(p)VG(x, q b (t, s))jS(x - q a (t, s)) d 3 x 

^ 0=1 Sa 6=1 56 

Thus, substituting the singular solution ansatz ( 12.71 ) and pairing with a 
smooth test function results in an expression in which (f) and its gradi- 
ent V0 appear linearly. 

Matching coefficients of (j) and V0 then yields a closed set of equations for 
the parameters p a (t, s) and q a (t, s), a = 1, 2, . . . , N, of the solution ansatz 
(ED, as 

(2.9) p a (t,s) = 0, 

5 N f 

(2.10) q a (M) = Kp)q~^2 Pb{sb)G(q a {t,s) , q b {t,s b ))ds b . 

a b—l 

Here the average density p is introduced in equation (12.81) . Thus, the density 
weights p a {t, s) = p a {s) are preserved, and the positions q a (t, s) in ( 12.71 ) fol- 
low the characteristics of the velocity ij = p(jj)V{G* p) along the Lagrangian 
trajectories given by x = q a (t, s). 

Remarks. 

• The result ( I2.9H2.10T ) holds in any number of dimensions. For example, 
it admits singular solutions of co-dimension one, supported on moving 
points in ID, along moving curves in 2D and on moving surfaces in 3D. 

• Holm and Putkaradze [3, [S] proved that the singular solutions ( 12.71 ) 
emerge spontaneously from smooth initial conditions and demonstrated 
their emergence in numerical simulations of equation ( 12.11 ) in one spatial 
dimension. In these ID simulations, the mass of each individual solution 
remained constant, as required by equation ( 12.91 ) . Moreover, when two 
singular solutions collided, they were found to add their weights pi and 

thereby "clumping" together. Eventually, all the singular solutions 
concentrated into a single "clumpon," whose weight (mass) equaled the 
total weight of the initial condition. 
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• The dynamics ( 12.9112. 101) for the singular solution ( 12.71) of the HP equa- 
tion ( 12.11) is a bit degenerate, because the weights p a (t, s) = p a (s) each 
turn out to be preserved. This occurs because substituting the singular 
solution into the HP equation for density produces only one term pro- 
portional to the test function, 0, which in turn yields trivial dynamics 
of the weights p a , a = 1, 2, . . . , N. The other terms are proportional to 
V0 and determine the dynamics of q a (t, s). 

• The general situation for the variational evolution of an arbitrary or- 
der parameter k G V (not just a density p G A n ) might be expected 
to possess multiple terms proportional to both and V0, and thereby 
produce nontrivial dynamics for both p a (t,s) G V and q a (t,s) G K 3 . 
Note, however, that the singular solution ( 12.71 ) of HP only existed, be- 
cause when paired with the test function the substitution of ( 12.71) 
into ( 12.11) produced no higher derivatives than V0. This was the key 
condition for possessing singular solutions of HP. 

2.3. Geometric order parameter equations. 

Cautionary note about signs: From now on, we will use the convention 
that the velocity of particles moving under the potential E is v = —pV5E/5p. 
This convention introduces a minus sign, so that the force will be defined as 
minus the gradient of the potential. Hence, the formulas for Darcy velocity 
in the previous section each acquire a "— " sign. The equations for continuous 
evolution in [71, [8] may be obtained by changing the sign of energy E — > —E 
in the final formulas. 

Geometry of order parameters. Order parameters for continua need not 
be densities. Instead, they may be amplitudes, or phases, or direction fields, 
or any other type of geometrical object, such as scalars, vectors, tensors, 
differential one-forms, two-forms, etc. Such geometrical objects are charac- 
terized by how they transform under smooth invertible maps with smooth 
inverses (diffeomorphisms, called diffeos for short). The corresponding in- 
finitesimal transformations of geometrical objects are defined as their Lie 
derivatives with respect to smooth vector fields. 

To create evolution equations for these quantities, we mimic the pattern of 
the gradient-flow equation ( 12.11) . Two features of this equation guide its gen- 
eralization: 
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First, we seek an evolution equation for a geometric quantity k based on the 
following conservation law in weak form, 



In other words, the geometric quantity k is Lie-dragged with velocity u[k], a 
vector field whose dependence on k must be defined as a suitable generaliza- 
tion of Darcy's velocity. 

Second, the vector field u which Lie-drags that object must be bilinear in the 
mobility \i and the variation of energy SE/5k. In the gradient-flow equation 
( 12.11 ), the Darcy vector field is — (//[k] V^) 11 . We seek the generalization of 
this vector field for an arbitrary geometrical order parameter. 

The previous derivation of the HP equation [3, [8] must be modified to admit 
an arbitrary geometric object k G V in any representation vector space V 
of the diffeomorphisms, rather than specializing to densities. Physically, this 
corresponds to identifying the order parameter by how it transforms under 
the diffeos and specifying its flux vector accordingly. This approach will allow 
us to address the formation of singularities in order parameters that are not 
densities. In particular, we compute the geometric variational flow for an 
order parameter k G V taking values in any vector space V, by following 
steps similar to those followed in [7], IE] for deriving HP. Namely, 



The first step invokes symmetry of the L 2 pairing and takes variations of 
the order parameter k G V by using the Lie derivative with respect to the 
vector field u{(j)). The second step introduces the diamond operation o. 
The diamond operation is defined in terms of the Lie derivative £ u action of 
a vector field u acting on variables a G V and b G V* that are dual under the 
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L 2 pairing by 

(2.14) {£ u a 5 fy = (u , a o b) = — (a , £ u b) = — (u , b o a) 

As before, the sharp (ft) operation raises vector indices in the key third step 
( 12.121) . which defines the velocity vector field u{4>) = — (0 o k)^ in terms of 
k G V, G V* and the diamond operation. This is the generalization of 
vector field rj in equation ( 12.21) from a density to an arbitrary vector quantity, 

(2.15) u(<f>) = -{cj) o k)K 

The last step used by ( 12.131 ) is the definition of the diamond operation to 
express the Lie derivative of the mobility £ u k with respect to the vector field 

5 V 

(2.16) « = (/iO^)». 

OK 

From the physical point of view, the vector field velocity defined by ( 12.161 ) 
provides the natural generalization we seek of the Darcy velocity fiS/SE/Sp. 
Thus we obtain the following geometric order parameter ( GOP) equa- 
tion, 

(2.17) _ = _£ (/iof) ^. 

When the order parameter is a density k = p G A", then diamond specializes 
to gradient, the Lie derivative becomes a divergence and one recovers the HP 
equation of [H IB]. Thus, the GOP equation ( 12.171 ) generalizes the concept of 
gradient flow of a density to "diamond flow" of any geometric quantity. The 
corresponding energy equation follows from ( 12.171 ) as 

dE /Ok 5E\ I 5E\ 

£ l ,. ~8E\»K, 



[2.18) = - M-L Iko 




~dt-\di'~5K~/ ~ \ *"("*E> ,ft " fo/ 

5_E 

5k 

In what follows, we shall consider several other forms of this equation for 
geometric order parameters in various vector spaces. The explicit forms of 
these GOP equations and their corresponding energetics need to be calculated 
using the definition of diamond and its properties. 

3. Properties of the diamond operation 



The diamond operation o is defined in ( 12.141 ) for Lie derivative £ rj acting 
on dual variables a G V and b G V* by 

(3.1) (6 o a , rj) = — (b , £ v a) =: — (b , arj) , 
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where Lie derivative with respect to right action of the diffeomorphisms on 
elements of V is also denoted by concatenation on the right. The diamond 
operator takes two dual quantities a and b and produces a quantity dual to 
a vector field, i.e., a ao b is a one- form density. For example, as we compute 
in more details below, if a is a scalar and b = 5d n x is an n form (where n is 
the dimension of the space), then 

a o b = (BVa) ■ dx ® d n x 

The o operation is also known as the "dual representation" of this right ac- 
tion of the Lie algebra of vector fields on the representation space V [T7] . 



When paired with a vector field 77, the diamond operation has the following 
three useful properties: 

(1) It is antisymmetric 

(b o a + a o b , 77) = . 

(2) It satisfies the product rule for Lie derivative 

( £^ (b o a) , 77} = ( (£^ b) o a + b o (£% a) , 77 }. 

(3) It is antisymmetric under integration by parts 

(db o a + b o da , ij) = . 

First, the o operation is antisymmetric, 

(3.2) (b o a , rf) = — (a o b , 77) , 
as obtained by using, 

(3.3) (b , £ v a) + (£ r] b , a) = , or, (b , a 77) + (b 77 , a) = , 
and the symmetry of the pairing 



Example. For example, if a = A • V and 77 = 77 • V are vector fields and 
b = B • dx <g) d n x is a one-form density in Euclidean coordinate notation, we 
find for the L 2 pairing in one dimension 

(a, b) = J a_\b = J (A-B)d n x 

In this case, integrating the product rule identity for contraction ( _l ) 

£ v (a J b) = (£ v a) _l b + a J £ v b 
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yields, for the density (3- form) £ v (a-ib) with homogeneous boundary con- 
ditions, 

j £ v {aAb) = J div((A-B)r7)d 3 £ = 

= J (£ v a) _l b + a _l £ v b 
= (£ v a , b) + (a , £ r] b) 

= (b , £ v a) + (a , £^6} = — (b o a , 77} — (a o b , 77} 

As expected, ( 13.21 ) follows for this case, re-affirming that the diamond oper- 
ation is skew-symmetric. In other notation, the intermediate equation may 
be rewritten as 

= (i^a , b) + (a , £ r) b) = ( ad ?? a , b) + (a , ad*6) , 

where £^6 = ad*6 for a one-form density b. In vector notation, 

£ 7] a = ad ?? a = [77 , a] = rj • VA — A • V?7 

and the intermediate equation is verified by integration by parts. 

Second, the o operation satisfies the product rule under the Lie derivative, 

(3.4) (£{ (b o a) , r)) = b)oa,n) + (bo (£t= a) , 77}. 

This property may be verified directly, as 

(£{: bo a, 77} + {bo £{= a, 77} = (6^77, a) - {bn£, a) 
= (a , b (ad^ 77)} = — (a o b , (ad^ 77)} 

(3.5) = - (adj(a o 6) , 77) = - (£ e (a o 6) , 77) = (^(6 o a) , 77) , 

upon using (&£, an) + (b^n , a) =0, implied by (13.31) . in the first step. 

Finally, the o operation is antisymmetric under integration by parts, 

(3.6) ( b o da , 77) = — (d& o a , 77) , 

as obtained from commutation of the two types of derivative (da) 77 = d{arj) 
and integration by parts, 

(3.7) (6, d(arj)) + (db, an) = 0. 

These three properties of o are useful in computing the explicit forms of the 
various geometric gradient flows for order parameters ( 12.171) . Of course, when 
the order parameter is a density undergoing a gradient flow, then one recovers 
HP from (pTTTf) . 
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4. DlSSIPATIVE BRACKET AND ITS PROPERTIES 

The formula for energy decay in ( 12.181) suggests the following bracket notation 
for the time derivative of a functional F\k], 



dF[hz] /Ok SF\ I n 5F' 



ST ~\ at 'to/ ~ \ ^f^'to 



5E\ / 5F 



5k J \ 5 k 
(4.1) =: {{E,F}}[k] 

Thus, the GOP equation (12.171 ) may be written in bracket form as 
(4-2) ^ = -£ {fl<:&K = {{E,K}}. 

The properties of the GOP brackets {{ E , F}} defined in equation ( 14.11 ) are 
determined by the diamond operation and the choice of the mobility /i[k]. 
For example, the GOP bracket in equation ( 14.11 ) is neither symmetric nor 
antisymmetric, for a general choice of mobility /j,[k]. However, any physical 
choice of mobility should produce strict dissipation of energy, i.e. 

(4.3) {{£,£}} <0 

One particular choice for the mobility satisfying this energy dissipation con- 
dition is /i[k] = kM[k], where M[k] > is a non-negative scalar functional 
of k. For this choice, one finds 

p,£}} = -lfW/ S o^ ( K og)'\<o 
4.1. Properties of the GOP Bracket: Leibnitz identity. 



Proposition 1 (Leibnitz property). The GOP bracket (\4-l\) satisfies the 
Leibnitz property [22] • That is, it satisfies 



{{EF, G}}[k] = F{{E, G}}[k]+E{{F, G}}[k] 
for any functional E, F and G of k. 

Proof. For arbitrary scalar functionals E and F of n and any smooth vector 
field 77, the Leibnitz property for the functional derivative and for the Lie 
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derivative together imply 



5k I \ 5k 

= E (^ S I- V ) +F (^ S I- V ) 



Choosing 77 = f k o j^j then proves the proposition that the bracket (14.11) is 
Leibnitz. 

□ 



4.2. Connection to Riemannian geometry. Following [12], we use the 



GOP bracket to introduce a metric tensor on the manifold connecting a "vec- 
tor" dfK and "co- vector" 5E/5k. That is, we express the evolution equation 
( 12.171 ) in the weak form as 



for an arbitrary element ip of the space dual to the k space, and where g K (- , •) 
is a symmetric positive definite function - metric tensor - defined on vectors 
from the dual space. (Here, we choose a + sign in front of metric tensor to be 
consistent with the choice of the sign for energy functional throughout this 
paper). 

Let us first notice that for any choice of mobility producing a symmetric 
bracket (in particular, /j,[k] = kM[k]) 

{{E,F}} = {{F,E}}, 

we may regard that symmetric bracket as defining an inner product between 
the functional derivatives, 

(4.5) {{ E,F }} =: gK (-,-) = U»-, (no 



5k, 
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Alternatively, ( 14.51 ) can be understood as a symmetric positive definite func- 
tion of two elements of dual space 0, ip: 

(4.6) 0«(<^) = ^o0, (Aco^y 

Notice that g((f), 0) > by ( Of ). 



Proposition 2 (Metric property of the GOP equation). For the choice of 
metric tensor ( [^.5|J^.6y ; £/ie GOP equation \2.11\ ) may be expressed as the 



metric relation H4-4\ ) 



Proof. For an arbitrary element ip in the dual space, 
/Ok 



= IftOlpi 

(4.7) = (,og, (K ov-)»)= ffB (f 




□ 



This approach harnesses the powerful machinery of Riemannian geometry to 
the mathematical analysis of the GOP equation ( 12.171 ). This opens a wealth of 
possibilities, but it also limits the analysis to mobilities fi for which the GOP 
bracket ( 14.11 ) is symmetric and positive definite, as in the modeling choice 
h[k] = kF[k]. Nevertheless, its strength seems to outweigh its limitations so 
we plan to apply this approach for modeling physical processes using GOP 
equations in our future work. 

5. Existence of singular solutions for the GOP equation (12.171 ) 

Let's choose free energy and mobility E[k], h[k] G H 2 and k G H 1 (which 
includes delta functions). For example, one may choose 5E/5k = G * k and 
mobility tensor fi\K] with k = H * k with Helmholtz kernels G and H with 
two different length scales, as done for the HP equation ( 12.11 ) from ( 12.11 ) in 

Big. 



The geometric order parameter equation ( 12.171 ) is then 



-^ = - £ 0iMo^ = - £ u[k}^, with u[k\= j 
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When paired with a smooth test function cf) dual to k, in L 2 , the solution k, 
satisfies, cf. equations ( 12.121 12.131) 




One may now substitute the singular solution Ansatz for GOP 
(5.2) K(x,t) = Pa(t,s)5(x - q a (t,s)) ds , 

a Js 

with k, p a £ V into the left hand side of ( 15.11 ). This substitution will produce 
only terms proportional to (j) and V0, after the appropriate integrations by 
parts. Fortunately, this precisely matches the corresponding terms on the 
right side for any k, since the Lie derivative for arbitrary tensor quantity (f) 
only contains terms proportional to (f) and its gradient. 



6. TWO FAMILIAR EXAMPLES 

We write the GOP equation ( 12.171) explicitly in two familiar cases. In each 
case we express the GOP equation of motion ( 12.171) in terms of the Leibnitz 
bracket ( 14.11) as 

(6-1) ^ = {{E , F }}[k] 

for any smooth functional F[k]. 

In three dimensional Euclidean coordinates, the Lie derivatives of a scalar / 
and a density D d^x are: 

(6.2) -£ v f = -v-Vf and - £ v (D d 3 x) = - V • (Dv) d 3 x . 

Having these two fomulas, we may compute from the definition of diamond 
(EH) that 

(6-3) ^[/]og = gv/i and Wjjj = "^ V S 
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6.1. Functions (bottom forms). For functions /, the relation ( 15.11 ) with 
diamond becomes, 

(6.4) (M <> g , (0 o /)») = - ( Vp[/]) , (0V/)«^ 

Consequently, the Leibnitz bracket ( 14.11 ) for functionals of a scalar function 
/is 

{{£\F }}[/] = - J g^V/ • V M [/]d 3 x 

As expected, this Leibnitz bracket is not symmetric for an arbitrary choice of 
mobility fi[f}. However, it produces proper energy dissipation {{ E , E }} < 
provided V/ • V/z[/] > 0. 

6.2. Densities (top forms). For densities, the diamond pairing is 

(6.5) (f»^(^*) , ) = (" V |- (flVtf ) 



and the Leibnitz bracket ( 14.11 ) is expressed as 



(6.6) {{E,F}}[D] = - J D M [D]V^| • V^d 3 x 



In this case 



J Tp 

^ = {{E,E }m 

Dfi[D] V 51 



d 3 x 



5D\ 

(6.7) < provided D/i[D] > 

Remark. These relations determine the dependence of the GOP energetics 
( 12.181 ) on the choice of mobility in other geometric cases case, simply by 
replacing (J) in fomula ( 14.51 ) above by the appropriate variational derivative 
of energy. 



7. GOP EQUATIONS OF MOTION 
The diamond relations above lead to the following equations of motion. 
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7.1. GOP equation for a scalar. First, for a scalar / 6 A 

df 

(7-2) = -£ { s_E vmi f 

(7-3) = -(^V/i[/]) tt -V/ 



Equation ( 17.11 ) can be rewritten in characteristic form as f(x,t) = const 
along 

dx 5E„ r „. 

This is an unusual characteristic equation, since the speed of characteristics 
depends on the nonlocal quantities 5E/5f and mobility fi. 

7.2. GOP equation for a density. For the density Dd 3 x G A 3 , we have 



7.4) = div^^plvg)') 



This recovers the HP equation (12.11) . (apart from the minus sign in front due 
to the revised definition of velocity, as noted in the beginning of section 2.3). 



8. Examples of singular solutions for the GOP equation (12.171 ) 
Substitution of the singular solution for GOP 

(8.1) K(x,t) = ^2 p a (t,s)S(x-q a (t,s))ds, 

a Js 

into both sides of ( 15.11) will produce only terms proportional to 4> and V0, 
after the appropriate integrations by parts. Matching corresponding terms 
yields the dynamics of p a (t, s) and q a (t, s). 
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For example, by GOP formulas ( I2.17U5.1U6.4I ). the scalar / satisfies, 
df\ 

' dt) 

■V/), (0V/i[/])t 



/og,^^!/]) 1 



= (/,div0(^V/z[/]) 



(8.2) 

One inserts 
(8.3) 



Sf 
■5E 



= (/, V0 • (^Vm[/]) + div (— V/4/]) 



/(M) = J] jpa(t,s)S{x-q a (t,s))ds 



with scalars p a and a = 1, 2, . . . , N into the left hand side of this equation 
and matches terms to find the equations for the parameters p a and q a of the 
singular scalar equation. These take the forms, 

'SE. 

Sf 



(8.4) 
(8.5) 



Pa(t,s) = p a (t, s) div ^— V/i[/] 
Pa(M) qa(M) 



' 5E 

Pa(t,s)(—Vfl{f] 



x=q a (t,s) 



One then chooses and //[/] so that the right hand sides make sense as 
functional relations. 



In Figure Q], we demonstrate the spatio-temporal evolution of a scalar given 
by (17.11) . We have taken 5E/5f = H * / where H is the inverse Helmholtz 
operator H(x) = e~^' a with length-scale, or filter width a = 1. The initial 
conditions consists of equally spaced (5-peaks in / with random strengths p a : 
/M) = £ a P a (0)£(x - q a (0)), -1/2 < p a (0) < 1/2. The evolution of the 
positions and amplitudes of the peaks is highly complex and shows strong 
sensitivity to the initial conditions. Thus, although the exact solution to 
the nonlocal PDE follows from the system of ordinary differential equations 
(I8.4U8.5I) . the solution of that system may still be highly complex and sensitive 
to the initial conditions of the p : s and g's. Fortunately, an exact analytical 
solution for this system of ODEs is available in the important case of evolution 
of a pair of ^-peaks. This solution imparts general understanding of the long- 
term behavior of (I8.4U8.5I) so we shall describe it in detail. 
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20 40 60 

Figure 1 . Numerical simulation demonstrates the evolution of an initial set of 8- functions 
in /. The vertical coordinate represents / = H * / ', which remains finite even when / 
forms 5-functions for our choice H(x) = e~' x '. This simulation of the initial value problem 
for equation (|7.ip uses generalized initial conditions f(x,0) = ^2 a p a (0)S(x — q a {0)) with 
random initial amplitudes —1/2 < p a (0) < 1/2. 

We assume initial conditions for the scalar / as 

f(x, 0) = p 5(x - q ) + sp 5(x - q ) . 

Here, the sign s = +1 specifies symmetric initial conditions, whereas s = — 1 
enforces antisymmetric initial conditions. We shall choose the sign s so the 
evolution due to ( I8.4U8.5I) preserves the symmetric position of the (5-function 
at x = ±q(t). If this choice is possible, then a solution of p(t),q(t) in terms 
of quadratures can be found. Not all energies E[f] and mobilities fi[f] allow 
symmetry preservation; however, many physically relevant cases do have this 
property. By reflection symmetry of the PDE, the even and odd parts of the 
solution are separately invariant. Consequently, we may assume a solution 
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for all t in the form, 
(8.6) f(x, t) = p(t) (5(x- q(t)) + s5(x + q(t)) 



Hence, the four equations (18 .408. 51) reduce by reflection symmetry to only 
two equations for p(i), q(t): 

(8.7) 



p(t) = p(t) div (— V//[/]) 



X=±q(t) 



6E 



V/4/] 



x=q a {t) 



Let us first notice that the derivative is taken at x = ±g, then 

x=±q(t) 



div 



x=±q(t) 



If we now denote 
(8.9) 



$>(p, q) 



■5E 
~5f 



V//[/] 



x=±q{t) 



then equations ( I8.TU8.8I ) can be written as 

d 



(8.10) 

We are now ready to prove the following 



Theorem 2 

Suppose ^(p,q) defined by (\8.9\) is homogeneous with ^(p,q) = p 1 i^{q), and 
s (symmetry parameter) is chosen so the solution retains its symmetry under 
the evolution. Then, a set of initial conditions f = po (5(x — qo) + s5(x + go)) 
exists, such thatq(t) in the solution ( Iff. 6]) collapses q — ► in finite time, when- 
ever ifj(q) is a continuous function bounded away from zero with \i^{q)\ > A > 
in some neighborhood of q = 0. 

Note. Many physical choices of energy and mobility admit the required ho- 
mogeneity of ^(p, q). For example [7], [8j selected 5E/5f = G*f,n = 1— H*f, 
where G and H are given functions. This choice of energy and mobility im- 
plies ^(p, q) = p 2 ip(q), so 7 = 2. Another example we shall employ here 

2 / 2 

is G(x) = e~ x ' a with a being a fixed parameter, and V pi = 1, since the 
explicit formulas are particularly simple. This choice yields 7 = 1. 
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Proof of Theorem 2 

First, let us notice that equations ( 18.101 ) integrate exactly (in terms of quadra- 



tures) if we assume that ^(p, q) = p 1 i>{q)- Indeed, ( 18.101 ) are equivalent to 

- . n dp n 1 d , , , 

which integrates exactly in terms of initial conditions p(0) = po, q(0) = qo as 



P = Po 



Notice that by the assumption \ip(q)\ > A > in some neighborhood of q = 
we can choose qo sufficiently close to so i/j{qo) 7^ 0. Then, the (/-equation of 
( KM gives 

(8.12) q = P^(q,y 2 ^(q? 1+l 

If, by assumption, ifj(q) < —A < if \q\ < 5 for some A > 0, 5 > 0, let us 
choose po > and < qo < 5. Then 

a < P ° A 27+1 
so q(t) goes to zero in time t* not exceeding 

(8 - 13) " 

The case ijj(q) > A > is treated analogously by choosing po < 0. 

We performed a simulation of the evolution of a scalar with initial conditions 
of the type ( 18.11 ). More precisely, we chose SE/df = G*f with G(x) = e~ x a , 

a 2 / 2 

and V/i = 1. This gave values 7 = l,ijj(q) = e 9 /a , so ^ satisfies the condi- 
tions of the theorem. We chose initial conditions / = po(5(x — qo) — 5(x + qo)), 
since in this case, the antisymmetric solutions preserve symmetry under evo- 
lution, with the choice po = — 1/8, qo = 5/3. 

The evolution of positions for ^-functions for this simulation is shown in 
Figure [21 The exact solution is shown with a solid line. The dashed line illus- 
trates the position according to the numerics. The exact solution collapses in 
finite time whereas the numerical solution shows exponentially slow approach 
of the £-peaks. We attribute the apparent discrepancy for large time to nu- 
merical dissipation. Indeed, introducing a term —ep mimicking numerical 
dissipation in the right-hand side of the p equation in ( 18.101 ) prevents solu- 
tion collapse, for any e > 0, however small e may be. Since every numerical 
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scheme must necessarily involve some numerical dissipation or distortion, we 
believe it may be difficult to construct a numerical scheme that shows exact 
collapse of the solutions. This numerical question will be addressed elsewhere 
in future studies. 




FIGURE 2 . Collapse of a solution with initial condition having two (5-functions with equal 
strength, opposite in sign /(x, 0) = po (S(x — qo) — S(x + qo)) with po = 1/8, qo = 2. Exact 
solution is shown with a solid line, results of numerics is given with a dashed line. 



The solution behavior may be investigated further by comparing the nu- 
merical and exact solutions for amplitude (18.41 ) in the evolution of a sin- 
gle #-peak. For this comparison, one may start with a single delta-function 
f(x,0) = PoS(x) and choose parameters ji = 1 — H * /, SE/Sf = G * / 
and G(x) = e— \x\, H(x) = \e~\ x \. In this case, equation ( 18.41) reduces to 
p = — 2j9 3 , whose solution is 

(8.14) 1 = ^ + 4* 

P 2 Po 
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t 

FIGURE 3. Evolution of the l/p(t) 2 versus time (circles). The theoretical prediction (|8.14[) 
is shown as a solid line. Note that (|8.14|) contains no fitting parameters. 

The equations for the singular-solution parameters p a and q a for the other 
quantities in equation ( IS. 3D may be found the same way, by substituting 
the solution ansatz ( 18 . 31 ) above into the formulas (I2.17U5.1U8.2I ) and matching 
terms. 

Not unexpectedly, the density case recovers previous results for the HP equa- 
tion. 

9. S0{2)- AND S0(3)-VALUED DENSITIES AND GYRONS 

9.1. Derivation of equations. An interesting situation for practical appli- 
cations occurs in the evolution of particles whose mutual interaction depends 
on their relative orientation. One example familiar from everyday life is the 
attraction of floating particles of non-circular shape, such as squares or stars. 
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In this situation, the attraction between any two particles arising through 
their mutual deformation of the surface depends on their relative orientation. 

Let us consider the motion of a set of particles whose energy E depends on 
their relative orientations. A reasonable physical assumption is that each 
particle carries its own measure of orientation. (This is the case with floating 
particles - stars and squares, say.) However, notice that the rotation angle 
and density are independent - one may find dense regions with little relative 
rotation and regions of small mass density but large relative rotation. This 
possibility leads us to study the evolution of a physical quantity k which 
carries mass density and orientation separately. Mathematically, this corre- 
sponds to a density (an element of A 3 ) that carries a "charge" taking values 
in the real numbers R for the mass and in the dual of the Lie algebra so(n)* 
with n = 2 or 3 for orientation in two or three dimensions. Such a quantity 
is written as n G (R <8> so(n)*) A 3 . These are (dual) Lie- algebra- valued 
"charge" densities. Let e a with a = 1,2, ...,n(n — l)/2 be a set of basis 
vectors for the Lie algebra so(n). These basis vectors satisfy 



and tf c are the structure constants of the Lie algebra so(n). For the case 
n = 3 one recovers the familiar result that tf c = e c ab is the completely anti- 
symmetric tensor density. Let us expand k in terms of the dual basis elements 



Here, the summation on the Lie algebra index a ranges from zero (for the 
mass density) to dim(so(n)) = n(n — l)/2 for the orientation degrees of 
freedom. Thus, a = corresponds to the real-valued mass density and 
a = 1, 2, ... , n(n — l)/2 denote the basis vectors of so{n). 

The Lie algebra so(3) has three basis vectors. In contrast, so(2) has only 
one basis vector, since only one number - angle of rotation - is sufficient to 
describe the rotation in two dimensions. Then, for example, in three dimen- 
sions we have k £ (R ® so (3)*) <S> A 3 . Correspondingly, the energy variation 
is given by 5E/5k G R (g) so(3) and is a function that takes real values for 
a = and takes values in the Lie algebra for a > 0. The distinction between 
so(3) and its dual so(3)* is purely notational for the present case, since these 
may both be identified with R 3 . However, we shall write the formulas below 



(9.1) 



[e b , e c ] = ijj c e<j (b, c, d > 0) 




(9.2) 



Kd n x = K n e a d n x . 
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in notation that would be valid for an arbitrary Lie algebra. 

The evolution equation for k 6 (10 so(3)*) A 3 in the space of densities 
taking values in the dual Lie algebra so (3)* is expressed as a pairing with a 
Lie algebra valued function G A°®so(3). The mobility fi G (R®so(3)*)cg>A 3 
lies in the same space as tz. Let us introduce the real-valued velocity vector 
r\ = — fi^D^ 6 , where the spatial covariant derivative operator D is expressed 
in components as, 

(9.3) Dj = dj + ad*^ , or (D j( f>) d = (dfo + t d ab A b )cf> a . 

Here, one sums over Lie algebra index b and the operation ( • )" in rj = 
— ftfrD^ 6 raises the spatial index. This vector field has components 

(9.4) rf = - Kb (D^) b :=- Kb (A/) * + Ka t- b A^ b . 

Thus, 7] is a real vector field in W 1 . This formula expresses Darcy law velocity 
for Lie- algebra- valued quantities by using the spatial derivative in covariant 
form, with Lie- algebra- valued connection A c e c whose j—th spatial compo- 
nent is Aj c e c , with j = 1, 2, . . . , n and (Aj C e c )$ = A? c e c . 

Likewise, we introduce the covariant time evolution operator as, 

(9.5) D* t K = dtK + ad^/c = (d t K a + t d ab A b K d )e a , 

where Aq = A^ea G so(3) is the temporal component of the Lie algebra valued 
connection. We also introduce the covariant spatial divergence, cf. equation 

as 

(9.6) Div*«?7 = Dj(Krf) = dj(Krf) + ad^Ar^ . 

Finally, we are ready to compute the covariant evolution equation. Applying 
the previous definitions and integrating by parts yields 

(D^,0) = (6k, ^) : =(- Div *^K^) = (M^D^) 

= ^M(-^DV), = (W («(/z[«]„D»||)) , <^ . 

Thus, we find the same equation we would have derived from thermodynam- 
ics, now written in terms of Lie- algebraic covariant derivatives, 

(9.7) DJac = Div* ^/x[/t] a D^jj . 
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As one might have expected, this is the covariant form of the HP equation 
(17.41 ) for the orientation "charge density" k = k^ taking values in the dual 
of the Lie algebra. This is the covariant evolution equation we sought for 
interaction energies that depend on both particle densities and orientations. 

The covariant form of Lie-algebra- valued scalar GOP equation. 

Comparing the patterns of equations ( 17.41 ) and ( 19.71 ) with that for the scalar 
case ( 17.11 ). we may immediately write the corresponding dynamics for a Lie- 
algebra-valued scalar / = f a e a in covariant form, as 

(9.8) D t f = -(^DM/]) tt -D/. 

The index notation for the spatial and temporal covariant derivatives of den- 
sities in the dual Lie algebra was explained, respectively, in equations ( 19.31 ) 
and ( 19.51 ). The scalar equation ( 19.81 ) may also be expressed in index notation 
by raising or lowering the Lie algebra indices appropriately. This observation 
also sets the pattern for generalizing other GOP flows to Lie-algebra values. 

Particular modeling choices for SO (3). The Lie-algebra- valued connec- 
tions Aq and A are not yet defined. One may identify a good candidate for 
Aq by examining the limit when the mobility vanishes, h[k] = 0. In this 
limit, when we choose the Lie algebra so (3), and make a particular choice of 
Aq, equation ( 19.71 ) reduces to the famous Landau-Lifshitz equation for spin 
waves, 

(9.9) D t *« = (dtK a - et b A b K d )e a = , 
or, in obvious vector notation, 

(9.10) dtK-A XK = 0. 

The Landau-Lifshitz equation now emerges upon making the choice Aq = 
6E/5k, so that 

5E 5E 

(9.11) d t K = — x k = k x Aft , when — = — Ak . 

OK OK 

It remains to choose the spatial connection A as our final modeling parameter. 
For the sake of simplicity, in the rest of the paper we shall make the choice, 

A • dx = dxx 1 £ so(3) 

with x m the order parameter group SO(3). This choice is a "pure gauge" 
connection whose spatial curvature is always zero. The quantity x £ SO(3) 
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is an auxiliary variable denned by relating its advective time derivative to the 
temporal part of the connection Aq as 

dx ^ ( 5E\I , , 5E 



OX ( Qej\* 

X '■= TT" + u ' ^X — AoX where u = [ k o —— and Aq = 
at v ok J 



5k 



From its definition, one finds an evolution equation for the spatial part of the 
connection form, 

(d X X' 1 )' = (A • dx)' = (d t A + u • VA + (Vu) T • A) • dx 

= oIAq + ad^ (A • dx) 
(9.12) = (djJ§ + ^Af)e a te 

where Aq = A^eb G so (3) is the temporal component of the Lie algebra valued 
connection. This formula closes the system. Given the choice Aq = 5E/8k 
and the relation Aq = XX f° r the auxiliary variable X-, this calculation of 
the evolution equation for A • dx = dxx 1 extends to any order parameter 
group. The resulting system is 




(9.13) D*k = Div* 

(9.14) d t A + u • VA + (Vu) T • A = V^ + ad Ao A 

where u = (^k o j^j and Aq = j^. In general, if the connection terms 

are small compared to spatial gradients, ( 19.141) becomes simply a nonlinear 
diffusion equation. Alternatively, if the connection terms are much larger 
than the spatial gradients, right-hand side of ( 19.141 ) becomes an integrable 
double-bracket equation considered in [21] )• The left hand side of ( 19.141) drives 
the evolution of k through the presence of time derivative and Landau-Lifshitz 
terms. 

9.2. Measure-valued solutions: gyrons. A special case that warrants 
further investigation arises when the Lie algebra in question is abelian, so 
that the structure constants vanish. As we shall show below, this case yields 
singular solutions. 

Let us multiply ( 19.71) by an arbitrary function G R <S> so(3). After the first 
integration by parts, we find 

(9-15) (°-^,A = UU\^ 



\dt' v ~\ \ n S 'dxiS K J 'dx 



J 
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By proceeding analogously to the calculation of HP clumpons in previous 
sections, we see that singular solutions of the form 

(9.16) K a = J p a (s, t)6(x - q(s, t))ds 

are admitted, in which p a takes values in the dual Lie algebra. Since the 
right hand side of ( 19.151 ) contains no terms proportional to a , we find that 
dtp a {s, t) = 0, for a = 0, 1, . . . , dim(so(n)). Equations for the j-th component 
of q which are expressed from the equation for component b are: 

^qj / (5E 

(9.17) p ^ =ps ^ MoDj — 

Assuming that p^ ^ 0, we get simple evolution equations for the strength p 
and coordinates qf 

(9-18) % = 9 4-UW 



dt dt V 5 k* 



(x=q) 



As we shall see, the singular solutions (gyrons) do appear in experiments. 
Thus, the singular solutions derived here are an essential part of the dynamics 
for particles with orientation. 

9.3. Comparison with experiments for the case of floating five-point 
stars, illustrate the power of our theory and connect it with the original 
problem of self-assembly of particles with orientation, we present a numerical 
simulation of the system (19.141 ) formulated for the particular case of a five- 
point floating stars. 

In order to formulate ( 19.141 ) explicitly, we must first find a suitable approx- 
imation for the energy as a functional of density and orientation. Denote 
the density at a point x as p(x) and orientation as a(x). If the energy be- 
tween particles due to density is given by the (half )-sum of binary interactions 
(which is a suitable approximation for dilute states), the 'density' part of the 
energy gives 

(9.19) Ep = \j P( X )P( X '^ G ( X _ x ') dxdx> 

For a given cr(cc), the orientation represents a vector with cartesian coor- 
dinates (coscr(cc), sincr(a?)). If we assume that a function H (a(x),a(bx')) 
describes interaction of orientations, then the energy due to orientations is 
given by 

(9.20) Ea = \j p ( x ^ p ( x ^ H M 3 ^' <7 ( aj/ )) dxdx ' 
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One possible approximation would be to assume that the interaction function 
K is simply proportional to the scalar product between <t(x) and a(x'). 
However, for the case of interacting stars or any objects with an m-fold 
symmetry, we must take into account that cr(x) and &(x') can differ by 
27r/m without affecting the interaction energy. Thus, the following form of 
the interaction energy is proposed: 



Here, Kir) is a scalar smoothing function, which must be determined from 
physical reasoning. An interesting case relevant for the evolution of the 
'pointy' objects is that the orientation part of the energy E a is of much 
longer range than the density part E p . The total energy is then simply 



where k = (p, a). 

In Figure [U we present a numerical simulation for the evolution in two di- 
mensions of a collection of five-point stars (m = 5). In this case, there are 
two unknown variables: the mass density p, which multiplies the basis vector 
of eg, and orientation a which is the coefficient of the basis vector e\. The 
initial conditions for the evolution are taken to be a gaussian distribution of 
initial density, possessing random initial orientations (Fig. [U top left). In 
the simulation, we observe the formation of large-scale structures which are 
predominantly either individual line segments or arrangements of three line 
segments meeting at a point to form a star-like structure (Fig. [H, top right). 
An experiment, conducted by P. D. Weidman, started with 4mm plastic float- 
ing stars distributed within a circular shape with random initial orientation 
(Fig. [U, bottom left). The stars are positioned on a net which is lowered 
slowly into the fluid. Care is taken to avoid residual water motion and con- 
vection. After 1 hour, the stars were found to assemble into line shapes or 
several-armed structures. Three-armed structures seem to be most common 
(Fig. [U, bottom right). While the position of the large-scale formations in 
the experiment is random and therefore impossible to predict analytically, the 
predominant shapes found in the experiment are reproduced in simulations 
using our model. 





(9.22) 



E[k] = E p + E t 
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H*p(t=0) H*p(t=0.1) 




10 10 




FIGURE 4. Top: Numerical simulation of density p = H*p where H(r) = e r is a Gaussian 

filter. White denotes the areas with high density whereas black marks the areas with small 

— i i 2 // 2 

density. Left: initial state for density, taken to be Gaussian:p(r, 0) = po.fee |r| ' p. In the 
simulation, the particles are assumed to have a five-fold symmetry. The initial conditions 
for orientation are random. Top right: final state. Bottom: an experiment reproducing this 
numerical solution (courtesy of P.D.Wcidman). Bottom left: initial conditions - five-point 
stars are positioned on a net and slowly lowered into the fluid. Initial orientations for the 
stars are random. Bottom right: final state after 1 hr. 
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